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Abstract. In continuation of [6] we discuss metrics of the form 

p 



(h,k) = f jh*i(Vol(f))3((Pi) f h,k)y6l(rB) 

on the space of immersions Imm(M, N) and on shape space £?; (M, N) = 
Imm(Af, AT)/ Diff(Af). Here (N,g) is a complete Riemannian manifold, M is a 
compact manifold, / : M — > N is an immersion, h and fc are tangent vectors 
to / in the space of immersions, f*g is the induced Riemannian metric on M, 
vol (/*<;) is the induced volume density on M, Vol(/) = f M vol(f*g), <3>; are 
positive real- valued functions, and (Pi) f are operators like some power of the 
Laplacian A-^ 3 . We derive the geodesic equations for these metrics and show 
that they are sometimes well-posed with the geodesic exponential mapping a 
local diffeomorphism. The new aspect here are the weights <£i(Vol(/)) which 
we use to construct scale invariant metrics and order metrics with positive 
geodesic distance. We treat several concrete special cases in detail. 



1. Introduction 

This paper is a continuation of the article [5] which discussed general Sobolev 
metrics on shape space: The special examples of metrics on spaces of shapes in l n 
studied there are invariant under the motion group, but not under scalings. The 
geodesic equation is well posed. 

In contrast, the article |4 j discusses weighted L 2 -metrics: They have easy hori- 
zontal bundles, and they can be made scale invariant by a judicious choice of the 
weight function. But we do not know that the geodesic equation is well posed. 

In this article we discuss in detail several interesting cases of weighted Sobolev 
metrics on shape space which are more explicit special cases of the general setup in 
[51 section 6]. In particular we are interested in scale invariant versions. For shapes 
of planar curves this was discussed in [T2]; a detailed study of a concrete example 
is in EH- 



1.1. The shape spaces in this work. Let (N,g) be a complete Riemannian 
manifold, and let M be a compact manifold (the template shape). A shape is a 
smoothly embedded surface in N which is diffeomorphic to M. The space of these 
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shapes will be denoted B e — B e (M, N) and viewed as the quotient (see [10] for 
more details) 

B e (M, N) = Emb(Af, N)/ Diff (A/) 

of the open subset Emb(M, N) C C°° (M, N) of smooth embeddings of M in N, 
modulo the group of smooth diffeomorphisms of M. It is natural to consider all 
possible immersions as well as embeddings, and thus introduce the larger space 
Bi = Bi (M, N) as the quotient of the space of smooth immersions by the group 
of diffeomorphisms of M (which is, however, no longer a manifold, but an orbifold 
with finite isotropy groups, see [TU]). 

Emb(M, N) ^ Emb(M, N) / Diff(M) B e (M, N) 



Imm(M, N) »- Imm(M, N) / Diff(M) Bj (M, N) 

More generally, a shape will be an element of the Cauchy completion (i.e., the 
metric completion for the geodesic distance) of Bi ( Af , TV) with respect to a suitably 
chosen Riemannian metric. This will allow for corners. In practice, discretization 
for numerical algorithms will hide the need to go to the Cauchy completion. 

Since Bi(M,N), as a quotient space or moduli space, is notoriously difficult to 
work on, all our computations will be done on the smooth infinite dimensional 
manifold lmm(M,N). In particular we shall use tangent vectors h to Imm(Af, N) 
with foot point / e Imm(M, N) which are smooth vector fields h : M —> TN 
along /. The tangent bundle TImm(M, N) is the smooth open submanifold of all 
h G C°°{M, TN) such that ir N o h is in the open subset Imm(M, N) C C°°(M, N). 

Any tangent vector h G T/Imm can be decomposed uniquely as h = Tf.h T + h , 
such that h^(x) is g-orthogonal to T x f.T x M at every point x G M. We call h 1 and 
h the tangential and normal parts of h, respectively. 

1.2. Sobolev metrics on shape space. As defined in [5J section 6], a Sobolev 
metric is a Diff (M)-invariant weak Riemannian metric on Imm(M, N) of the form 

Gf(h,k)= / g(Pfh,k)vol(g) for / G Imm, h.keTflmm. 
J M 

where g = f*~g is the induced Riemannian metric on M, vol(<?) = vol(f*g) is the 
induced volume density on M, and where P is a smooth section of the bundle 
i(TImm; Tlmm) over Imm which is Diff ( M )-invariant such that at every / G Imm 
the operator 

Pf : Tyimm — » Tjlmm 

is an elliptic pseudo differential operator of order 2p that is symmetric and positive 
with respect to the 7f°-metric on Imm, 

H° f (h,k)= [ g(h,k)vo\(g). 
Jm 

(1) If moreover the smooth adjoint 

Adj(VF) G T(L 2 (Tlmm; Tlmm)) 
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of VP exists in the following sense: 

/ g((V m P)h,k)vol(g)= [ g(m, Adj(VP)(/i, fc))vol(ff), 

JM JM 

then the geodesic equation for the Sobolev metric exists, see [5J 6.2-6.5]. 

(2) In [51 theorem 6.6] we have proved that the geodesic equation for a Sobolev 
metric is well posed in the sense that we have local existence and uniqueness for 
geodesies in the C°°-sense in such a way, that the Riemannian exponential mapping 
is a local diffcomorphism, if P satisfies the following further hypothesis: 
P, VP and Adj(VP) 1 - are smooth sections of the bundles 

L(TImm;TImm) _L 2 (TImm; Tlmm) L 2 (TImm; TLmm) 



Imm Imm Imm, 

respectively. Viewed locally in trivializations of these bundles, 

p f h, (VP) f (h,k), (AdKVPMM)^ 

are pseudo-differential operators of order 2p in h, k separately. As mappings in the 
foot point / they are non-linear, and it is assumed that they are a composition of 
operators of the following type: 

(a) Local operators of order I < 2p, i.e., nonlinear differential operators 

A(f)(x) = A(x, V l f(x), V'- 1 /^), . . . , Vf(x), f(x)), 

where V is a fixed covariant derivative independent of /, and 

(b) Linear pseudo-differential operators of degrees Z j , 
such that the total (top) order of the composition is < 2p. 

For numerical experiments concerning Sobolev metrics of order one see the arti- 
cles is □]. 

1.3. Almost local metrics. In the paper [1] we studied metrics of the following 
type on the shape space S,(M, M. n ) (where dim(M) — n — 1): 

G*(h, k)= f 3(Vol(/), Tt(Lf)) g(h, k) vo\(f*g), f € Imm, h, k € P/Imm, 

JM 

where now g is the Euclidean metric on R", vo\(f*g) is the volume density on 
M induced by the the pullback metric g = f*g on M, Vol(/) = J M vo\(f*g) is 
the total (n — l)-volume of /(M), L* is the Weingarten mapping of the immersed 
hypersurface, Tr(L^) is the mean curvature, and where $ : K 2 — >• M.>o is a smooth 
function. Note that Tr(L^) is a nonlinear differential operator of / and thus local. 
By calling such a metric almost local we meant that the only non-local expression 
in the metric is the total volume Vol(/). The investigation in 4 was aimed at 
curvature computations and numerical examples. Doing numerics for almost local 
metrics is facilitated by the fact that the horizontal bundle for almost local metric 
is very simple. In the paper [5] this investigation was extended to metrics of the 
form 

Gf(h,k) = f H>(Tr(Lf),det(Lf))g(h,k)Yol(tg) 

JM 

where now $ is a positive function of the mean curvature and the Gauss curvature. 
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1.4. The metrics in this paper: weighted Sobolev metrics. Here we combine 
the settings of the papers [6], [4], and [5] in the following way: The ambient space 
(N,~g) is again a complete Riemannian manifold, but not necessarily R". The 
metrics are of the form 

(1) G f (h,k)= [ 5><(Vol(/)) g{{P^ f h,k)v6L{rg), 

Jm s :=o 

where the Pi are smooth sections of the bundle i(TImm; TImm) over Imm which 
are Diff (M)-invariant such that at every / 6 Imm each operator 

(Pi) f : Tf Imm — > Tjlmm 

is a pseudo differential operator that is symmetric and positive semidefinite with 
respect to the i?°-metric on Imm, where the <E>i are positive smooth functions, 
and where J2i ^i(Vol(/)). (Pi)f is elliptic and positive self adjoint. The operator 
governing the metric is thus Pf := ^i(Vol(/)).(Pj) / . We may again call the 
corresponding metric G p an almost local metric if all (Pi)f(h) are local operators 
in h and in /, i.e., differential operators according to the linear ([13], [14]) and 
non-linear ([171) Peetre theorem, and if the only non-local ingredients are functions 
ofVol(/). 

In section [2] we discuss general metrics of the form (1), deduce the geodesic 
equation (12.31 if the adjoints of VPi all exist) and remark that it is well-posed 
(under the conditions listed in II. 2\i , where we allow the assumptions (|1.2I 2) in a 
more general form: The linear operator vol(/*g) n- J M vol(/*7?) = Vol(/) is not 
a pseudo-differential operator in the strict sense since £(£) is never an admissible 
symbol. But the operator is very easy to handle and thus [6] theorem 6.6] continues 
to hold. In 12.41 we discuss the conserved quantities arising from symmetry groups. 
Finally we discuss the horizontal bundle and the horizontal geodesic equation which 
corresponds to the geodesic equation on shape space Bi(M,N). 

Theorem 12.61 discusses conditions which imply that geodesic distance on shape 
space is positive. We recall the conditions for positivity given in [6] and add two 
conditions for operators of order zero. This is a nontrivial question, since geodesic 
distance on shape space vanishes for the L 2 -metric (where P — Id) as shown in [11] 
and [10) . and also on diffeomorphism groups (see |10j ) and on the Virasoro group 
(see [5]). We even have vanishing geodesic distance of diffeomorphism groups for 
Sobolev metrics of order < p < \, see [3]. 

Section [3] is devoted to almost local metrics of order in h, where the operator 
is given by P f = £\ $i(Vol(/))(P i ) / = $i(Vol(/)).^(/). Id where *,-(/) is a 
smooth positive function depending smoothly and equivariantly on /; for example, 
\&(/) = || Tr 5 (5)|| 2 can be the g = /*g-norm squared of the vector valued mean 
curvature. For these metrics the horizontal bundle at / is just the bundle of vector 
fields along / which are g-orthogonal to /. In particular, we discuss the case of con- 
formal metrics where Pf = $(Vol(/)). Id, and of curvature weighted metrics where 
P = (1 + A\\ Tr 9 (S)\\^) Id. We derive the geodesic equations for all these metrics. 
Note that the general well-posedness theorem [UJ theorem 6.6] is not applicable for 
these metrics since the order in h is not high enough. 
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Section [4] is devoted to metrics corresponding to operators 

P f = f]Vol(/)^- 1 3((A0 l / l ,fc)vol(r5), 

i=0 

where A? is the Bochner-Laplace operator associated to / and g, acting on vector 
fields on M along /, see [SI 3.11 and 5.9]. The powers here are chosen in such a 
way that for N = W 1 the metrics are invariant under rescalings. We derive the 
geodesic equations on the space of immersions and on shape space. If p > 1 then 
the geodesic equation is well posed. 

The final section [5] is an appendix which recalls the induced covariant derivative 
V 9 on any bundle over the space Imm(M, N) for the convenience of the reader, 
and we differentiate the vector valued mean curvature Tr g (S^) with respect to the 
immersion /. 

2. Weighted Sobolev type metrics 

We want to study metrics that are induced by operators of the form 
P f ^^^(VolCf)) (Pi)/ : T/Imm(M, N) ^ T f Imm(M,N). 

Assumption. For each i £ I let : K + i-> R + be a smooth function and Pj be a 
smooth section of the bundle of linear maps i(TImm; TImm) over Imm. In addition 
we assume that Pf is an elliptic pseudo differential operator that is symmetric and 
positive with respect to the _ff°-metric on Imm, i.e.: 

Hj(P f h, k) = Hj(h, P f k) and Hj(P f h, h) > 0, for h, k e T / Imm(M, AO- 
Note that an elliptic symmetric operator is self-adjoint by [TCI theorem 26.2]. 
Then the operator P induces a metric on the set of immersions, namely 

Gf(h,k) := / s(P,M)vol( 5 ) = 5>i(Vol) / g{(Pi) f h,k)vol(g) , 
Jm ieI Jm 

for / G Imm, and h, k E Tjlmm. The metric G p is positive definite since P is 
assumed to be positive with respect to the H°-metric. 

2.1. Invariance of P under reparametrizations. 

Assumption. It will be assumed that all Pi are invariant under the action of the 
reparametrization group Diff (M) acting on Imm(M, N), i.e. 

p. = (r<pyp. f or a ii ^ e DifT(M). 

For any / £ Imm and if G Diff(Af) this means 

[Pi)f = {T f r*)- 1 o(P i ) foi/)0 T f r*. 

Applied to h G Tjlmm this means 

(Pi) /(h) oip= (Pi)f 0l p(ho ip). 

Note that Vol is invariant under the action of the diffcomorphism group and there- 
fore all (Vol) are invariant as well. This together with the invariance of Pj implies 
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that the induced metric G p is invariant under the action of Diff (M). Therefore G p 
induces a unique metric on Bi as explained in [HJ section 2]. 

2.2. The adjoint of VP^ . Following [Bj we introduce the adjoint operator of VP^ to 
express the metric gradient H which is part of the geodesic equation. Hf arises from 
the metric Gf by differentiating it with respect to its foot point / G Imm. Since the 
metric is defined via the operator P = JT (Vol). Pi, one also needs to differentiate 
(Pi)/ with respect to its foot point /. As for the metric, this is accomplished by 
the covariant derivate. For P, G r(L(TImm; TImm)) and m G TImm one has 

V m P t G T(L (TImm; TImm)), WP t G r(i(T 2 Imm; TImm)) . 

See [51 section 2] for more details. 

Assumption. For all Pi there exists a smooth adjoint 

Adj(VP 4 ) G T(T 2 (TImm; TImm)) 
of VPi in the following sense: 

/ g((V m Pi)h,k)vo\(g)= [ 5(m,Adj(VP,)(M))vol( 5 ). 

JM JM 

The existence of the adjoint needs to be checked in each specific example, usually 
by partial integration. 

Lemma ([HI)- If the adjoint of VP^ exists, then its tangential part is determined 
by the invariance of Pi with respect to reparametrizations: 

Adj(V Pi)(h,k) T = (g(VP^,fc)-5(V^,P 4 fc)) J 

= grad 9 g(Pih, k) - (g(P t h, VA;) + g(Vh, P ( fc)) J 

for f G Imm, h,k G I/Imm. 

2.3. The geodesic equation on the manifold of immersions. As explained 
in the introduction the geodesic equation can be expressed in terms of the metric 
gradients H and K (see [6j section 4.4] for more details on this topic). These 
gradients and thus also the geodesic equation will be computed in this section. 

Theorem. If the adjoint Adj(VP;) exists for all i £ I, then also H and K exist 
and the geodesic equation is given by: 

V d Jt = \H f {f u ft)-K f (f t J t ) 

= 1 P ~ X (T, ^(Vol^AdjCVP.XA, f t ) x - J^9{Pifu ft>ol(g) Tr 9 (S) 

- 2Tf.g(P i f t ,Vft) i -9(P i fuft).Tr 3 (sj)) 
-^ -1 (l>i(Vol)(V A Pi)/t-y; / .9(/ t ,Tr 9 (5))vol( 3 ))^(Vol)PJ f 
+ Tr* (g(yf t ,Tf)).Pf t ) . 
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Remark. Note that we can apply [6J theorem 6.6] . Therefore the geodesic equation 
is well-posed assuming that the operator P satisfies all conditions listed in [6l section 
6.6], which are reviewed in (|1.2I 2). 

Remark. As shown in [B] the geodesic equation can be equivalently written in terms 
of the momentum p — G(/t, .) as: 

P = G(f t ,-) 

V 9t P=^G/(ff (/*,/*),•), 

For the metric G p the momentum takes the form 

p = Pft® vol( 5 ) = $(Vol)P/ t (g> vol(g) : K -> Tlrnm 

since all other parts of the metric (namely the integral and~g) are constant and can 
be neglected (cf. [BJ section 6. 5] ). Then the geodesic equation reads as: 

p = 5^* < (Vol)P i / t ®vc%) 



= \ ( E ( Vo1 ) ( Ad i ( v ^) (/* . /*) 

*i(Vol) 



- 2Tf.g(P i f t , Vf t f - g(Pif u ft).^(S) ) ) ® volfo) 



Proof of Theorem \2.3i The governing equations for the metric gradients H,K £ 
r(L 2 (TImm; TImm)) are 

(V m G)(h, k) = G(K(h, m),k)= G(m, H(h, k)) , 
where h, k,m are vector helds on Imm (cf. [BJ section 4.3]). We calculate: 



(V m G )(h, k) = D(/, m ) / g(Ph,k)vo\(g) 



(1) 



M 

g{P(V m h),k)vo\{g)- [ g(Ph,\7 m k)vol(g) 

M JM 

D(f, m )9(Ph, k)vol(g) + / g(Ph, k)D {f , m) vo\(g) 

M JM 

g(P(V m h),k)vo\(g)- [ g(Ph,V m k)vol(g) 

M JM 

g(V m (Ph),k)vo\(g)+ [ g(Ph,V m k)vo\{g) 

M JM 



g(Ph,k)D (ftm) vol(g) 



M 

g(P{V m h),k)vol{g)- [ g(Ph,V m k)vol{g) 

M JM 

= [ g{(y m P)h,k)vol(g)+ [ g{Ph,k)D u . m) vol(g) . 

JM JM 

Since P is an elliptic, self-adjoint, and positive operator, it is invertible on the space 
of smooth sections. (See the beginning of the proof of [BJ theorem 6.6] for a detailed 
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argument.) Plugging in the variational formulas from [5] for the volume form and 
the volume, one immediately gets the K- gradient. 

K f (h,m) = p- 1 (V$ i (Vol)(V m P i )ft-^ I 9(m,Tr s (S'))vol( ff ))^(Vol)P i /i 

i Jm 

+ Tr s (g{Vm,Tf)).Ph) . 
This formula is equal to the formula presented in [6] section 6.3]. 

To calculate the ff-gradient, one rewrites equation ([T]) using the definition of the 
adjoint of P (c.f. 0| section 6.3]): 

(V m G p )(h,k)= f Vv m ($ J (Vol)POM)vol( 5 )+ f g(Ph,k)D {f , m) vol(g) 
Jm i Jm 

= V^Vol) / g(m,Adi(VPi)(h,k))vol(g) 
Jm 

+ y)$5(Vol)D (/im jVol / g(Pih,k)vol(g) 
Jm 

+ V$ 4 (Vol) f g{P l h,k)D {Lm) vo\{g). 
Jm 

Now the second and third summand are treated further using the variational for- 
mulas for the volume density and volume from section [5] For the second term this 
yields: 

V ^(Vol^^Vol f g(Pih, fc)vol(g) = 
Jm 

= -y>S(Vol)/ g(m,Tr 9 (S))vol{g)) [ g(PA k)vo\(g) . 
Jm Jm 



For the third term we calculate: 



X>i(Vol) / g(P i h,k)D {f>m) vol(g) = 
Jm 

= V $,(Vol) / g(Pih, k) Tr* (g(Vm, Tf))vo\(g) 
Jm 

= V$,(Vol) f g(Pih,k) Tt 9 (Vg(m,Tf) - g(m,VT f))vol(g) 

= 5>,(Vol) f g(P t h,k)(-V*g( m ,Tf)-g(m,Tr°(S)))vo\(g) 
Jm v ' 

= -V$ 2 (Vol) / 5 ?(V5(P^,fe),3(m,T/))vol(p) 
Jm 

-V^(Vol) f g(PAk)g(m,Tr°(S))vo\(g) 
Jm 

= V ^(Vol) / g(m, -Tf. grad 9 g(Pih, k) - g(P t h, k) Tr 9 (S))vo%) 
Jm 
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$5 (Vol) 



g(Pih,k)vol(g)Tr°(S) 



Collecting terms one gets that 
Gf(H f (h,k),m) = (V m G p )(h,k) 

= W^Vol) / s(m,Adj(VPi)(M) - , [V / 

- T/. grad 9 g(Pih, k) - g(P, I, . k ) Tr« (S))^o^ !1 ) 

Thus the iJ-gradient is given by 
ff / (/ l ,fc) = p- 1 ('^$ i (Vol)(Adj(VP i )(/ l ,fc) 



*i(Voi) y M 



/ g(PA fc)vol(g) Tr 9 (S) - T/. grad 9 s(PA *) - S(P^, *)• Tr 9 (S) 



The highest order term grad 9 g(Pih, k) cancels out when taking into account the 
formula for the tangential part of the adjoint from section 12.21 



H f {h, k) = p- 1 f MVo\) ( Adj(vp)(/>, ky 

$'fVol) f 

' 1 g(ph,k)vol(g)^(S) 



**CVol) J M ' 

Tf.(g(Pih,\7k) +g(Vh,P i kj) t -g(Pih,k).Tr 3 (S))] . □ 



2.4. Conserved Quantities. The metric G is invariant under the action of the 
reparametrization group on M. According to [4] the momentum mapping for this 
group action is constant along any geodesic in Imm(M, N).: 

VX g X(M) : G(Tf.X, f t ) rep. mom. 
or G(/ t T , •) g r(T*M ® M vol(M)) rep. mom. 

If (N,g) has a non-trivial isometry group, this gives rise to further conserved quan- 
tities with values in the Lie-algebra of Killing vector fields. For example, take a flat 
ambient space N = R™ . Then the almost local metrics are in addition invariant un- 
der the action of the Euclidean motion group R™ x SO(n). This yields the following 
conserved quantities: 



p 

5^$i(Vol(f))Pi/ t vol(s) hn. mom. 

i=l 

VX g so(n) : / V $i( Vol(f) )s(P/ t) X/) vol( 5 ) ang. mom. 
or / V $i( Vol(f) )(/ A P/ t ) vol(p) g A 2 »" = *o(n)* ang. mom. 



/ 
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2.5. The horizontal bundle and the geodesic equation on Bi. According to 
[6] the geodesic equation on Bi is equivalent to the horizontal geodesic equation 
on Imm assuming that in every equivalence class of curves there exists a horizon- 
tal curve. Therefore we need to study the horizontal bundle of the metric G p . 
Multiplication with an everywhere positive function does not change the notion of 
horizontality. Therefore we can repeat the analysis of [6j sections 6.8-6.9] to obtain: 

Lemma. The decomposition of h G Tf Imm into its vertical and horizontal compo- 
nents is given by 

h"° r =(Pj)- 1 ((P f h) T ), 

h hor = h- Tf.h ver = h- Tf.iPjy 1 {{P f h) T ). 

Furthermore for any smooth path f in Imm(M, N) there exists a smooth path if in 
Diff(M) with (p(0, . ) = Mm depending smoothly on f such that the path f given 
by f(t,x) — f(t,tp(t,x)) is horizontal: 

G p (dJ,Tf.TM)=Q. 

Thus any path in shape space can be lifted to a horizontal path of immersions. 

See [6l sections 6.8-6.9] for a proof of the lemma. Note that the above lemma es- 
tablishes a one-to-one correspondence between horizontal curves in Imm and curves 
in shape space. Thus the geodesic equation on shape space is given by: 

Theorem. The geodesic equation for the G p -metric on shape space is equivalent 
to the set of equations for a path of immersions f : 

p = Pf t ®vol(g), P/ t = (P/ t ) x , 



( v ^) hor = ^(E $ '( Vol )( Ad j(VA)(/t,/t) 

^ i 

^i(Vol) 



(I)(V()]i/M ^/ t ,/ t )vol( 5 )Tr fl ( 5 ) 

-5(PJ f ,/ t )-Tr 9 (S)) ) ®volfo) 



The proof is a direct consequence of theorem 12.31 and of lemma 12.51 

Theorem 2.6 (Geodesic distance). The metric G p induces non-vanishing geo- 
desic distance on B e if at least one of the following conditions is satisfied for 
Ci, C2, C3, A > and for all h G TImm: 



(i) 




>Ci\\h\ 


(2) 




> c 2 voi 


(3) 




>C 3 \\h\ 



Proof. For the proof of (1) we refer to [H section 7]. Assuming that (2) holds we 
calculate for a horizontal path /: 



T Imm 

Len G p 



(/)>C 2 \ WVol(/) f g(ft,ft)voKg)dt>C 2 [ f l.Jg(f t J t )vo\(g)dt 

JO V JM Jo JM v 
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= C2 / v °l(/(-7 -)*g) = C2 (area swept out by /) . 

J[Q,l]xM 

Since area swept out is a positive distance on B e this shows the statement for (2). 
For condition (3) we calculate: 



Len^TV) > C 3 J^ ^Jjl + A\\ T^(SW)g(f tl f t )vo\(g)dt 

I \ [ g(ftJt)vo\(g)dt>C 3 f 1 -L= f l.Jg(f t J t )vo\(g)dt 
h V J M Jo V Vol Jm v 



>C 3 
> C 3 



mWVol(/(t,.)) 



(area swept out by /) 



(3) will be proven by showing that y/Vo\(f(t, .)) is Lipschitz continuous. Using the 
formula from section 15.21 for the variation of the volume we get 



d t Vol(/) = - / g(f u Tr 9 (S))vo\(g) < 



A I 

{ I l 2 vol(.g) 

' M 



g(f u T^(S))vo\(g) 



A I 



g(f t ,Tr°(S)) 2 vol(g) 



M 

<VWT){ I \\Ml\\Tr°(S)\\lvol(g) 



A I 



4> 



<yVol(7) / ||/ t |Hvol(g) < ^M[f)JG f {f u f t ) 



M C2 



ICt. 



Thus 



v y > 2/voi(7) " 2V^V njt'jv 

By integration we get 



VWhj- VvoiCM- / 5 tv ^oT(7)dt 



< 



1 -±= y jG f (f t ,ft) = ^-L=Len« m (/). □ 



3. Almost local metrics of order zero 



In this chapter we will study metrics that are induced by operators of order zero, 

i.e. 

p = £ $(Vo\(f)).(p f Uh) = ]T $cvoi(/)).*i(/).fc, 

where ^ : Imm(M, AT) t— > M>o is any smooth, positive function depending on the 
immersion /, that is equivariant with respect to the action of the diffeomorphism 
group Diff(M). Note that ^ might also be a differential operator acting on the 
immersion /. So order zero means order zero in h, fc, but possibly higher order in 
the foot point /. The most prominent examples of such operators, and also the 
examples we will consider in this section are 

P(f) :=$(Vol)vI/(||Tr2(S)|| 2 ). 
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Here Tr 9 (5) denotes the vector valued mean curvature. 

For the case of planar curves, metrics of this form have been introduced in [12] . 
Recently they have been generalized to hypersurfaces in n-space, see [4] . The most 
general case of surfaces in a possibly curved ambient space has been studied in the 
PhD thesis of Martin Bauer [7]. 

The great advantage of this class of metrics is the particularly simple form of the 
horizontal bundle, namely: 

Lemma (Horizontal bundle). For an almost local metric of order zero, the hor- 
izontal bundle at an immersion f equals the set of sections of the normal bundle 
along f. 

Proof. By definition, a tangent vector ft, to / G Imm(M, N) is horizontal if and 
only if it is G-perpendicular to the Diff (M)-orbits. This is the case if and only if 
~g(h(x),T x f.X x ) = at every point x £ M. □ 

This observation is of particular interest both for numerical and theoretical rea- 
sons. The simple form of the horizontal bundle allows to solve the boundary value 
problem for geodesies on shape space numerically by minimizing the horizontal en- 
ergy (c.f. [4]). And it allows for the computation of the sectional curvature, as 
done in [12] for the case of planar curves and in [4] for the case of hypersurfaces 
in n-space. The main drawback of this class of metrics is that it is unclear if the 
geodesic equation is well posed. (The operator P does not satisfy the conditions of 
[3 theorem 6.6].) 

In the following we want to study two examples of almost local metrics in more 
detail. 



3.1. Conformal metrics. The class of conformal metrics correspond to metrics 
G p with 



For the case of planar curves these metrics have been treated in [18] , pjj] , [20] , [15] 
and for the case of of hypersurfaces in n-space they have been studied in [4] . In the 
case of planar curves [TS] provides very interesting estimates on geodesic distance 
induced by metrics with <f>(Vol) = Vol and e Vo1 . 

According to theorem 12.51 we can simply read off the geodesic equation both on 
the manifold of immersions and on shape space. It is given by: 

Theorem. The geodesic equation on the manifold of immersions for the class of 
conformal metrics is given by: 



If one passes to the quotient space P>i (M, N) the geodesic equation reduces to: 



P = $(Vol).Id . 



P 



$(Vol)/ t ® vol(g) 




p=*(Vol)/t®vol(ff), 



ft = ft e Nor(/), 
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(Va t p) hol ' = ~(V(Vol) / g(f t J t )vo\(g)Tr*(S) 
z v Jm 

+ <P(Vo\).g(f t J t ).Tr°(S)) ®vol( ff ). 

Remark. Note that the geodesic equation for conformal metrics does not change 
when passing from a flat to a curved ambient space. 

3.2. Curvature weighted metrics: The G^-metric. In our setting the curva- 
ture weighted metric G A corresponds to the operator 

P=(l + A\\Tr g (S)\\ 2 ) . 

Note that P is a differential operator of order two applied to the foot point /, but 
an operator of order zero in h and k. Metrics weighted by curvature have been 
studied in [11], [12], @], [5]. 

According to section 12.51 we need to calculate the adjoint of P to obtain the 
geodesic equation. 

Lemma. For the operator P = 1 + A\\ Tr 9 (5)|| 2 the adjoint is given by: 
Adj(VP<)(fc> k) = 4A Tr (g^Sg-^S, Tr 9 (S))g(h, k)) 

- 2A(A( Tr 9 (S).g(^ fc))) X - 2A§(h, k) Tr 9 (iF(T/, TV^Tff 
+ ATf.g(h 1 k).gr^\\^{S)\\l 

Proof. Using the variational formula for the mean curvature (see [5]) we calculate: 

f g{m,Adj(VPi)(h,k))vol(g)= [ g((V m Pi)h,k)vol(g) 
Jm Jm 

= a[ (D ftm \\Tr3(S)\\ 2 )g(h,k)vo\(g) 
Jm 

= 2A f g(TT\S),D Lm Tr\S))g(h 1 k) V o\(g) 
Jm 

= 2A f g(2TT(g(m ± 7 g- 1 .S.g- 1 ).S),T^(S)).g(h,k)yo\(g) 
Jm v 7 

-2A f (g® 3o °)(A(m ± ),TrS(5)).3(/ i ,fc)vol( 3 ) 
Jm 

+ 2Aj^g(Ti-V(RV(m ± ,Tf)Tf),Tr°(S))g(h,k)vo\(g) 
+ 2Aj \d ||Tr s (5)||| (m T ).g(h, k)vo\(g) 
= 4A^Tr(3(m x , ff - 1 .S. 5 - 1 ).3(5,lV s (5)))3(/i,fc)vol( ff ) 
- 2Aj (g ® g° ) (^(m ± ),'It a (S))^h, k)vol(g) 
+ 2Aj TV (g(BF(m x ,Tf)Tf,T^{S)))g(h,k)vol(g) 
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+ A J <?(grad s ||TrS(^)||| ! m T ).^,fc)vol(.g) 
AA g (m 1 - , Tr (<T 1 .S.g~ 1 .3(5, Tr« (5)) *) J ) v, ,1 ( .„ i 
(g ® ff0 °) (m x , A ( Tr* (5).^, fc)) ) volfo) 



-2A 



M 

-2 f ATr° (g(R°(Tf, Tr 9 (S))Tf, m x ))g(h, k)vol(g) 
+ Aj g(Tf. gradf ||T^(5)||i, T/.m T ) .g(h, k)vol(g) 
AA g (m 1 , Tr («T 1 .S.g- 1 .g (S, Tr 9 (S)) .g(h, fc)) ) vol( 5 ) 
-2A / yfm- 1 , A(Tr 9 (S*).g(/i,fc)))vol( 5 ) 
-2 f g(g(h,k)ATr 9 (^(T/, Tr«(5))T/), m ± )vol( 5 ) 
+ ^^g(™,r/.3(/ l ,fc).grad 9 ||Tr3(5)|||)vol(. g ) □ 



Using the above lemma and theorem 12.51 yields the geodesic equation for the 
G^-metrics both on the manifold of immersions and on shape space: 

Theorem. The geodesic equation for the G A -metric on Imm(.M, jV) is given by 

p=(l+A||Tr 9 (S)|| 2 )/ t ®vo%) 
V 9tP = i^ATr^-^-^S.Tr^S^M)) - a(a( Tr*(S).s(M)) ) ^ 

- fc) Tr« ( Ji*(T/, Tr^ (5)) Tf) ^ 
-{l + A\\^{S)\\ 2 )gUuft).^{S) 

- 2(1 + A|| Tr 9 (S)\\ 2 )Tf.g(f t , V0) ® vol( 5 ) 



TTie geodesic equation for the G A -metric on shape space is equivalent to the set of 
equations for a path of immersions f : 

p= (1 + A|| Tr 9 (S)\\ 2 )f t ® vol(a), / t = G Nor(/), 
(Va tP ) hor = ^2ATr( ff - 1 5.g- 1 5(^T^(5))^,fc)) - a(a( Tr*(S).g(M)) ) ^ 

- ^(/j, fc) Trs (Tf, Tr 9 (^))T/) ± 

- 1(1 + A|| Tr^(5)|| 2 ) g(f t , ft), ® vol(g) 
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4. Scale invariant Sobolev metrics 

In this section we will study scale invariant Sobolev type metrics of order p. This 
class of metrics is induced by an operator of the form 

p p 
P = J2 Vol" » ^A* = ^^(VoijPi, 

i=0 j=Q 

with 

$i(Vol) = Vor^ 11 - 1 and P = A\ 
Recall that m denotes the dimension of the base space, i.e. dim (A/) = m. 

Remark. There are other choices of scale invariant metrics. For example, one 
could use integrals over curvature terms as weights. A scale invariant metric of 
order zero that induces positive geodesic distance is induced by the operator: 

P = Vor-- 1 + Vor 1 ||Tr 9 (,S , )|| 2 . 

However, in this paper we decided to focus on volume weighted Sobolev metrics since 
they seem to be the most natural examples. 

In order to calculate the geodesic equation we need to calculate the adjoint 
operators. 

Lemma. The adjoint of VPj defined in section \2.2\ for the operator Pj = A 1 is 
given by: 

i-l 

AdKVA^XM) = 2^Tr(.g- 1 ScT 1 3(VA l - z - 1 /i,VA / fc)) 
1=0 

i-l 

+ (V'stVA- 1 - 1 /!, A l k)) Tr 9 (5) 

1=0 

i-l 

+ Tl ' 9 (i^CA*- 1-1 ^ VA l k)Tf) 

1=0 

i-l 

- Tr 9 (FfVA'-'- 1 /!, A l k)Tf) 

1=0 

+ Tf.[grad g g(A i h,k) - (g(A% Vk) + g(Vh, A**))*" . 

Proof. The proof of this lemma is similar to the calculation of the adjoint for P = 
1 + AA P in [61 lemma 8.2]. □ 

Using this we can calculate the geodesic equation for scale invariant Sobolev type 
metrics. 

Theorem. The geodesic equation for the weighted Sobolev-type metric of order p 
on shape space is equivalent to the set of equations for a path of immersions f: 
p 

P = Pft = J2 Vol" » "^Vt ® vol( 5 ), Pf t = ft € Nor(/), 

i=0 
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i=0 ^ 1=0 

i-l 

+ E (V^tVA'-'" 1 /*, A'/ t )) Tr 9 (S) 
;=o 
i-i 

+ £ Tr 9 (i?9(A 4 -'- V t , VA l f t )Tf) 

1=0 
i-l 

- £ Tr 9 (^(VA'- 1 - 1 /^ 1 /^/) 

-4/ A / (Ai/t ' /t)v ° 1(5)Tr9( ^ 

-g(A7t,/ t )-Tr 9 (5))) ® vol( 5 ) 

For p > 1 i/ie operator P satisfies all conditions of [6l theorem 6.5] and therefore 
the geodesic equation is well-posed. 

5. Appendix: The covariant derivative over immersions and some 

variational formulas 

In this appendix we will first recall a concept that has been introduced in [6], 
namely a covariant derivative on the manifold of immersions that is induced by the 
metric g on the ambient space N. This covariant derivative is used to calculate 
the metric gradients and to express the geodesic equation both on the manifold of 
immersions and on shape space. In the second part we will present some variational 
formulas, mainly without proof. For the missing proofs see for example [4], [6], [5], 

i, m- 

5.1. Covariant derivative on immersions. This section is taken from [6]. Let 
V 9 be the covariant derivative on N with respect to the metric g. This covariant 
derivative induces a covariant derivative over immersions as follows. Let Q be a 
smooth manifold. Then one identifies 

h G C°° (Q, TImm(M, N)) and X e X(Q) 

with 

h A 6 C*°°(Q x M, TN) and {X, M ) € X(Q x M). 

Using the covariant derivative 

Vf x , 0M) h A eC°°(QxM,TN), 

one can define 

V x h = (v 9 X0m) / 1 a ) V e C°°(Q,TImm(M,N)). 

This covariant derivative is torsion- free (see (6j). It respects the metric ~g but in 
general does not respect G. 
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It is helpful to point out some special cases of how this construction can be 
used. The case Q = E will be important to formulate the geodesic equation. The 
expression that will be of interest in the formulation of the geodesic equation is 
Va t /t, which is well-defined when / : R — > Imm is a path of immersions and 
ft : ffi — > TImm is its velocity. 

Another case of interest is Q = Imm. Let h,k,m €E 5E(Imm). Then the covariant 
derivative V m /i is well-defined and tensorial in m. Requiring V m to respect the 
grading of the spaces of multilinear maps, to act as a derivation on products and to 
commute with compositions of multilinear maps, one obtains a covariant derivative 
V m acting on all mappings into the natural bundles of multilinear mappings over 
Imm. In particular, V m P and V m G are well-defined for 

P e r(L(TImm; TImm)) , G G r (Z^ ym (TImm; K)) 

by the usual formulas 

(v m p)(h) = v m (p(h)) - p(y m h), 

(W m G)(h, k) = W m (G(h, k)) - G{V m h, k) - G(h, V m k). 

5.2. Variational formulas. Recall that many operators like 

g = f*g, S = S f , vol(<7), V = V 9 , A = A 9 , ... 

implicitly depend on the immersion /. In this section we will present their derivative 
with respect to /, which we call the first variation. The following lemma contains 
a compilation of results that can be found for example in [5] , [5] . 

Result. Let f be an immersion and ft € Xyimm a tangent vector to f . We have: 

(1) Let F : hxaa(M,N) r(TJM) be a smooth mapping, that is equivariant 
with respect to pullbacks by diffeomorphisms of M . Then the tangential 
variation of F is its Lie- derivative: 

(2) The differential of the pullback metric f*g is given by 

D {fJt) g = 2 Symg(Vf t) Tf) = -2g(f t x ,S) + 2 Sym V(/ 4 T ) b 
= -2g(ft,S)+C f T 9 . 

(3) The differential of the inverse of the pullback metric is given by 

D U,h)9^ = D U , ft) {rg)- X = 2g(f t L ,g- 1 Sg- 1 ) + C f7 (g- 1 ). 

(4) The differential of the volume density is given by 

D(f,f t )Vol(g)=Tr° (3(V/ f ,T/))vol( 5 ) - ( div 9 (/ t T ) - g(f^ Tr* (S))) vol( 5 ). 

(5) The differential of the total Volume is given by 

D(f Jt ) Vol(/) = D [fJt) I volte) = - / g{ft^{S)) vol( 5 ). 
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(6) For A e r(L(TImm; TImm)), h € T/Imm one has 

(V„A)(„) -tt(rtlW*-*») - V (v,o„„„» )+ ^ 1Wm „„ a ,)''' 
+ V(R'(f,,Tf)h) - Tr« (/?(/,, T/)Vfc 

For the proof of the above formulas in a similar notation we refer the reader to 
[51 13] • The following lemma is concerned with the first variation of the vector valued 
mean curvature: 

Lemma. The differential of the vector valued mean curvature 
Imm -> r(Nor(/)), / i-> Tr s (5) 

is given by 

D( fJt) Tr°(S) = Tv(2g(f t ± , g - 1 .S.g- 1 ).S) ~ A(f t x ) 

+ Trs XR9(f t x , Tf)Tf) + C f T Tr»(S). 

Proof. By formula 1 of the above result, the formula for the tangential variation 
follows from the equivariance of the vector valued mean curvature form with respect 
to pullbacks by diffeomorphisms. For the variation in normal direction we calculate: 

d t Ti 9 {S) = d t Tr(g-\S) = Tr ((d t g^).S) + Tr (g-^dtS) 

Using the formula for the variation of the inverse metric it remains to calculate 
the derivative (in normal direction) of the vector valued shape form. By definition 
S(X, Y) = V x {Tf.Y) - Tf.VxY- Therefore we have 

d t S f (X,Y) = V dt V x {Tf.Y) - V 9t Tf.V x Y 

= VxVvTfA + B?{Tfd u TfX)TfY - V VxY Tf.d t 
= V x V Y f t - V Vx y/t + RP{ft,T fX)T fY 
= V 2 f t +R°(f t ,Tf)Tf, 

where we used the following two formulas for interchanging covariant derivatives 
(see [6j): 

v (8t) o M )T/.(o R ,y)-v (QB) y ) r/.(ft J o Jl f) = o 

and 

v ( a t ,ov)V (Q ,,y ) ft-v (0B , y) v (8t ,o M) /i = fl: y (T/.(ft,Oj lf ),r/.(0R,y))/i. 

Thus we get 

d t Tr 9 (5) = d t Trig-KS) = Tr ((dtg-^.S) + Tr (g-'.dtS) 
= Tr {2g(f t ± ,g- 1 Sg- 1 ).S) + Tr (g-\V 2 f t x ) 

+ Tr (g^I&tff, Tf)Tf) + % Tr*(S) 
= Tr{2g(f t ± ,g- 1 .S.g- 1 ).S)-A(ft ± ) 

+ Tr (g- l BF{ft,Tf)Tf) + C f r Tr 9 (S) □ 
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